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We examine a model of M-component quantum rotors coupled by Gaussian- 
distributed random, infinite-range exchange interactions. A complete solution is 
obtained at M = oo in the spin-glass and quantum-disordered phases. The quantum 
phase transition separating them is found to possess logarithmic violations of scaling, 
with no further modifications to the leading critical behavior at any order in 1/M; 
this suggests that the critical properties of the transverse-field Ising model (believed 
to be identical to the M — ► 1 limit) are the same as those of the M = oo quantum 
rotors. 
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Extensive attention has been lavished in the last decade on the problem of classical spin 
glasses and their finite temperature (T) phase transition to a spin glass ordered phase [|TJ. 
In contrast, there has been relatively little work on quantum spin glasses [[2]|$],|J, especially 
on their T = quantum phase transition to a quantum disordered (or 'spin-fluid') state. 
In particular, there is no model for which the critical properties of this quantum phase 
transition have been obtained. On the experimental side, there has been a renewed interest 
in a number of spin systems which are in the vicinity of a T = phase transition from a 
spin-glass to a spin- fluid state HHIllllil • these include the dipolar, transverse- field Ising 
magnet LiHo x Y\- x F^ ||, the lightly-doped cuprates @,0,|J, and various layered transition- 
metal/rare-earth oxides [|]. 

In this paper we examine a quantum spin glass which allows us to examine more carefully 
the nature of the quantum spin-glass to spin-fluid phase transition and determine the spec- 
trum of excitations in the spin-fluid phase. We consider M-component quantum rotors with 
Gaussian-distributed random, infinite-range exchange interactions. A complete solution of 
this model will be obtained at M = oo in both the spin-glass and spin-fluid phases and at 
the critical point separating them. We also examine the nature of the 1/M corrections at 
T = in the spin- fluid phase and at the critical point: we find that the form of the leading 
critical behavior and the low-frequency spectral weight remains unmodified to all orders 
in 1/M from the M = oo result. Thus the results of this paper could have been derived 
without any reference to the 1/M expansion, by simply resumming Feynman graphs which 
are dominant at low frequency - these graphs happen to be identical to those selected by 
the M = oo theory. 

The quantum rotors should not be confused with true quantum Heisenberg spins present 
in any isotropic antiferromagnet; the different components of the rotor variables all commute 
with each other, unlike the quantum spins. As a consequence, the path-integral written in 
the rotor variables has an action which contains no Berry phases and is purely real. The 
properties of random quantum spin models are quite different from those of the quantum 
rotors considered here, and will be discussed elsewhere [fLl]] . Apart from its theoretical 
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simplicity, the main utility of the rotor model is that the M = 1 limit of the path integral 
is expected to be in the same universality class as the Ising model in a transverse field. 
The absence of any 1/M corrections noted above, suggests that the critical-behavior of the 
infinite-range, transverse-field Ising model is identical to that of the M = oo limit solved 
in this paper. This is also consistent with a recent analysis of this Ising model by Huse 



and Miller fll2| : their results for the critical point are essentially identical to those obtained 
below in the M = oo model. 

We will study the following ensemble of Hamiltonians 

H=^-YL 2 + ^!=YJ il h i -n j h\ = \ (1) 

where i,j extend over N sites, n ifl are the M components of a unit-length rotor on 
site i, the L ifJLV (/jl < v, fj,, v — 1 . . . M) are the M{M — l)/2 components of the angular- 
momentum generator Lj in rotor space, and the J^- are mutually uncorrelated exchange 
constants selected with probability P(J{j) ~ exp(— J^/(2 J 2 )). The are mutually com- 
muting variables and the quantum dynamics is defined by the commutation relations: 

The Lin U satisfy the commutation relations of angular momenta in M dimensions. As g — > 0, 
the model reduces to the classical, infinite-range, M-component, Heisenberg spin glass which 
was analyzed earlier by de Almeida et.al. [p70f| . 

The formulation of the N —> oo limit of H can be obtained by a straightforward general- 
ization of the analyses in Refs 0,|3],|J . We use the path-integral formulation of the partition 
function, introduce n replicas, and average over the ensemble of the J^. The iV — >• oo 
limit yields a saddle-point which describes the quantum mechanics of n replicas of a single 
rotor. Assuming the saddle-point is O(M) invariant (this is true in both the spin-fluid and 
spin-glass phases) we obtain the single-site path-integral 

Z = j Vh a {r)5 (h a2 (r) - l) exp f -— ^ dr{d T h a ) 2 + ^ drdT'Q ah {T - r> a (r) • h b {r 

(3) 



and the self-consistency condition 

Q ab (r - r') = (n a (r) ■ h\r')) Zo (4) 

Here a,b = l...n are replica indices, r, r' are Matsubara times, and (3 = 1/T. The 
Edwards- Anderson order parameter |TJ for the spin-glass phase is 

Qea = Q aa (r -> oo). (5) 

Moreover, Q ab , a ^ b, is r-independent and non-zero only in the spin-glass phase [|]]. 

An exact evaluation of Zq is clearly not possible. We present below the results of a 
systematic 1/M expansion on Z . 

M = oo theory: Imposing the constraint by a Lagrange-multiplier A, the M = oo limit 
of Eqns (|3]||) reduces to the constraint Q aa {r = 0) = 1 and 

Q{iuJ n ) = g (ul + A - gJ 2 Q{iu n )) 1 (6) 

where Q{iuj n ) is the Fourier transform of <3(t) at the Matsubara frequencies, and the r.h.s 
is a matrix inverse in replica space. 

1. Paramagnetic phase: For large g, or large T, we expect a paramagnetic phase (the 
quantum-disordered phase is the T = paramagnetic state) in which case Q ab will be 
replica diagonal H,|J . A closed- form solution can be obtained from @ for the spectral 
weight = Im(Q aa O + i0+)): 

»/ \ / J(a; 2 -A + 2^)(A + 2J( ? -^)] 1/2 
X (w) = sgn(w) (7) 

for A — 2Jg < uo 2 < A + 2Jg and x" = otherwise. It is clear that a physically sensible 
solution requires A > 2Jg where A is determined by the constraint equation n a2 = 1, or 

\ — x>)coth03w/2) = 1 (8) 
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It is evident from (|^) that the M = oo paramagnet has a gap of (A - 2Jg) 1 ' 2 towards spin- 
wave excitations. We expect 1/M fluctuations to fill in this gap at any finite T; the gap in 



the T = spin-fluid phase is however robust towards such corrections. The paramagnetic- 
spin glass phase boundary is determined by setting A = 2 Jg and solving (|8|) for a line in the 
g — T plane: the results of this calculation are shown in Fig [I]. The quantum transition near 
T = occurs at g = 9tt 2 J/16 — 3T 2 /J + ■ ■ •, and the classical transition near g = occurs 
at T = J — g/12 + • • •; this latter result agrees with that of Ref. (TH . 



2. Spin-glass phase: We now expect only Q ab {iu n = 0) to acquire off-diagonal compo- 
nents the finite-frequency Q(iu n ) remains diagonal. We therefore parametrize 

Q aa {iu n ) = Q a r a eg (iu n ) + /3q EA 5 Wn>0 (9) 

where Q™ g (iui n ) can be obtained immediately from the solution of (0) and continues to 
have spectral weight Xregif^) which obeys (0) with a value of A to be determined below. 
We parametrize the off-diagonal components of Q (iu n = 0) by an arbitrary hierarchical 
matrix |HJ specified by a monotonic function f3q(x) on the interval < x < 1. Using 
the expressions for the inverse of an hierarchical matrix in Ref. |TJ|], the self-consistency 
equation @ can be transformed into two integral equations for q{x) and q E A- Simple 
algebraic manipulations then yield the satisfactory [0] result 

Qea = q(l) (10) 

Repeated differentiation of the integral equations showed that dq/dx = 0; q(x) can therefore 
only be a piecewise constant function. We chose q(x) = q\ for < x < u and q(x) = q E A 
for u < x < 1, whence the integral equations specified q\ = and q E A] u wa s however 
left undetermined It was then necessary to evaluate the free energy and demand 

stationarity with respect to u. The final result was quite simple: we found u = implying 
that q(x) = q E A for all x and that the replica-symmetric solution is optimal. This agrees with 
the classical limit at g — which was found in Ref. [TIJ to possess a stable replica symmetric 



solution at M = oo; we also undertook a stability analysis, similar to that in Ref. [^0J, for 
the quantum-rotor model and found only non-negative eigenvalues in the fluctuations about 
the replica-symmetric state. Our final results for the spin-glass phase were: A = 2Jg with 
Xregi^) & Yen by (0) being gapless over the entire phase, Q ab {iu n = 0) = /3q E A for a ^ b and 
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q E A = l- — x^Mcoth(/?u;/2) (11) 
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3. Quantum critical region: We now examine the region near the quantum phase transi- 
tion at g = g c = 9tt 2 J/16, T = 0. Scaling (see e.g Ref. ||) predicts that the spin-glass 
paramagnetic boundary obeys T ~ \5g\ zu (here Sg = g — g c )- From the equation for the 
phase-boundary at small T above, we deduce zv = 1/2. The order-parameter q E A must 
vanish as q E A ~ I^Pi from ([11]) this yields (3=1. Further the T = spin-gap, A, in the 
quantum-disordered phase should vanish as A ~ (Sg) zu . Using A = (A — 2Jg) 1 ^ 2 and (||), we 
find however that A ~ (5g /\og(l/ 5 g)) 1 1 2 . Thus there is a surprising logarithmic violation 
of naive scaling - the log divergence is a consequence of the square-root threshold in the 
spectral weight (0). For u and 5g small, but uj/Sg arbitrary, the entire T = 0, local dynamic 
susceptibility obeys a scaling form: 

x "(uj,T = 0) = dsgnHM^ i^j (12) 

where the frequency scale A 9 obeys A g = C2(Sg) zv / log 1 / 2 (1/5^) for small Sg, the exponent 
\i = — 1 + (3/{zv) = 1 ||, Ci,c 2 are non-universal constants, and $ 9 is a universal function 
given by 

[ (1 - l/x 2 ) 1 / 2 for \x\ > 1 
*,(*) = (13) 
otherwise 

We will argue below that the results for zv, (3, n and $ 5 are in fact exact to all orders in 1/M; 
only the non- universal constants c\,c<i get modified by higher order corrections. A related 
analysis can be performed at the critical coupling g = g c but at finite temperature 0. For uo 
and T small, but with uo/T arbitrary, the local dynamic susceptibility now obeys the scaling 
form: 

XV> 9 = 9c) = cisgn(^) \uj\^ t (^-^j (14) 
where the universal function $^ is 
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;i-47r7(3* 2 ))V2 for |*| > 2n/VS 
$ r (x) = { (15) 

otherwise 

and the frequency scale A T = k B T / log 1//2 (l/T) at low T, with no non-universal prefactor. 
Note again the presence of logarithmic violations of naive scaling; the frequency- scale for 
the dynamic susceptibility is however still set completely by the absolute temperature to 
leading-log accuracy. The presence of a gap in $t is clearly an artifact of the large M 
limit [H], as the T = state is gapless at g = g c ; we expect 1/M corrections to modify $r 
by filling in the gap. 

1/M expansion: We now examine corrections to the above mean field theory at T = 
in the quantum-disordered phase and at the quantum-critical point, g = g c . We will not 
examine such corrections in the spin-glass phase where the structure is considerably more 
complicated due to the expected appearance of replica symmetry breaking. Our main result 
will be that neither the critical exponents nor the form of the low frequency spectral weights 
are modified by the 1/M corrections. We begin by absorbing all higher-order corrections 
into a self-energy, E, in the n propagator, which modifies @ to 

Q(ioj n ) =g(uj 2 n + \- gJ 2 Q(iu n ) + ^{iu n )/My l . (16) 

The function E(r) is itself a non-linear functional of Q(r), obtainable by a 1/N expansion of 
Zq. Let us consider first the critical point g = g c and use the M = oo result Q aa {iuj n ) ~ \uj n \ 
at low frequencies. The leading term in E satisfies Im(E(u; + i0 + )) ~ uj 5 , Re(E(u; + i0 + )) ~ 
ai+ci2U) 2 , at small oo; the suppression at low-frequencies in Im(E) arises from restriction in the 
phase space to three spin-wave decay. On the imaginary frequency axis, this implies that the 
leading non-analytic term in E(ia; n ) is ~ \oo n \ 5 . Now consider the self-consistency (|16|). The 



analytic terms in E lead to apparently innocuous frequency and mass renormalizations, while 
the non-analytic terms vanish so rapidly that they don't modify the assumed low-frequency 
form Q aa (iuj n ) ~ \ui n \; our initial assumption is therefore self-consistent. Terms higher-order 
in 1/M have even weaker non-analytic contributions. Thus there are no modifications to the 
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critical properties, order- by-order to all orders in 1/M. Similar considerations also apply to 
the low-frequency spectrum in the spin-fluid phase, where again the M = oo form survives. 

This behavior can be better understood using a classical statistical mechanics point 
of view, in which the system is viewed as a classical one dimensional spin system with 
a long range interaction Q aa (r); having solved the model we can then require the self- 
consistency (gj). Our results above imply that the critical point of the quantum phase 
transition corresponds to a spin system with Q{r) ~ 1/r 2 for large r (1/t 2 is the Fourier- 
transform of |o;|). We may consider a lattice discretization of r, and also replace the fixed 
length spins by M-component soft spins S with a Landau- Ginzburg potential local in time. 
Thus we are led to a model with action whose continuum limit is 



S = - J drdr' Q(t - t')S(t) ■ S(r') 

(17) 



+ 1 I dr 
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-(d T S(T)) 2 + rS 2 (r)+u(S 2 ) 2 
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where g,r,u are constants. This classical spin system, with Q(r) ~ l/r 1+CT , was studied 
many years ago These authors found a high temperature paramagnetic phase with 



power-law spin correlations, and a transition to a low-temperature ordered state if a < 1, or 
if M = 1, a = 1. In the high temperature phase they found (§(t) ■ S(r')^ ~ l/r 1+cr which is 
also the result obtained from the leading term in the high-temperature expansion (expansion 
in powers of Q) . Throughout the high-temperature phase the spin-spin correlation exponent 
is unmodified by higher order terms. As Q(t) and (S(t) ■ S(r')^ have the same asymptotic 
decay, it is evident that the self-consistency (|j) can be satisfied for any value of a. The 
result that a = 1 corresponds to the quantum phase transition, can be traced to the u; 2 in 
(§) or flT7|) which is generically present as the leading analytic u n dependence. It thus has 
nothing to do with the critical point of the one dimensional system; the quantum-critical 
point corresponds to a point in the high-temperature phase of the classical spin model. 
With the choice a — 1, the other critical properties then follow; the logarithmic violation of 
scaling comes in this model from summing tadpole diagrams in the 5* 4 interaction. These 
arguments are valid for all M including M = 1 (the transverse Ising case). 
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This paper has presented a soluble model with infinite-range interactions which displays 
a quantum phase transition from a spin-glass to a spin fluid phase. Further theoretical work 
on the extension of these results to finite-range interactions is clearly required. Our result 
for }("{uj) in (|7|) has qualitative similarities to the experimental results in the transverse-field 
Ising magnet LiHo x Y\- x F± |5[], and the outlook for an eventual complete understanding of 
this system is bright. 

We thank A. Georges, D. Huse, J. Miller, and A. P. Young for useful discussions. This 
research was supported by NSF Grants Nos. DMR 8857228, DMR 9157484 and by the A.P. 
Sloan Foundation. 
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FIGURES 

FIG. 1. Phase diagram of H in the T-g plane at M = oo. The line g = corresponds to the 



classical model of Ref [10]. The quantum-disordered phase is the paramagnet at T = 0. Regions 
in which the spin fluctuations are primarily thermal/quantum are noted. 
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